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Abstract 

We calculate the nucleoli parameters in isospin asymmetric nuclear matter using 
^vQ . the QCD sum rules. The nucleon self-energies S^, and S* are expressed in terms of the 

K^ I in-medium values of QCD condensates. The simple approximate expressions for the 

self-energies are obtained in terms of these condensates. Relation between successive 



(N 

^^ . inclusion of the condensates and the meson-exchange picture of the nucleon interaction 

^O I with medium is analyzed. The values of the self-energies and of the symmetry energy 

Z^ ■ agree with those obtained by the methods of nuclear physics. 

^ I. INTRODUCTION 

-*— > ■ 

^ ' In this paper we investigate the vector and scalar self-energies of nucleons in nuclear 

C ■ matter composed by the neutrons and protons, distributed with densities p„ and pp. We 
calculate the dependence on the total density 



X 



P = Pp + Pn (1) 

and on the asymmetry parameter 

Pp + Pn 



P = ^^— ^ (2) 



(another conventional presentation of the asymmetry parameter is /3 = ^r^ = 1 — ^ with 
A^ and Z standing for the total number of neutrons and protons, while A = N + Z). We 
present also the equations for the single-particle potential energies. The results are obtained 
by QCD sum rules approach. 

The QCD sum rules were invented by Shifman et al. |1^ to express the hadron parameters 
through the vacuum expectation values of QCD operators. Being initially used for the 
mesons, the method was expanded by loffe j2 El to the description of the baryons. The 
approach succeeded in describing the static characteristics as well as some of the dynamical 
characteristics of the hadrons in vacuum — see, e.g., the reviews jHE]- 



The main idea is to consider the correlation function 

Mq) = g,7"ng(g2)+/n^(g2), (3) 

(/ stands for the unit 4x4 matrix) describing the propagation of the system with the 
quantum numbers of the hadron, in the different regions of values of the momentum g^, 
where certain information on its behavior is available. The asymptotic freedom of QCD 
enables to present Uliq"^) (j = q,I) at g^ -^ — oo as the power series of g~^ and QCD 
coupling as- On the other hand, the imaginary part of UqIq'^) at g^ > can be described 
in terms of the observable hadrons. This prompts to consider the dispersion relation for the 
function Ili{q^) [T] 

at g^ —>■ — oo. The coefficients of the expansion of the left-hand side (Ihs) of the functions 
nQ(g^) in powers of g~^ are the expectation values of the local operators constructed of 
the quark and gluon fields, which are called "condensates". Such presentation, known as 
the operator product expansion (OPE) provides the perturbative expansion of the short- 
distance effects, while the nonperturbative physics is contained in the condensates. The usual 
treatment of the right-hand side (rhs) of Eq. (JH) consists in "pole + continuum" presentation, 
in which the lowest laying pole is singled out while the higher states are approximated by 
the continuum described by the pure perturbative OPE contribution. Thus, Eq. (JH) ties 
the values of QCD condensates with the characteristics of the lowest hadronic state. Such 
interpretation requires that the contribution of the pole to the rhs of Eq. (j^ exceeds the 
contribution of the continuum. 

The OPE of the Ihs of Eq. Q becomes increasingly true, while the value of |g^| in- 
creases. On the other hand, the "pole + continuum" model becomes more accurate while 
|g^| decreases. The important assumption is that the two presentations are close in certain 
intermediate region of the values of g^. To improve the overlap of the QCD and the phe- 
nomenological descriptions, one usually applies certain mathematical anzats, i.e., the Borel 
transform. The Borel transformed dispersion relations Q are known as QCD sum rules 

una- 

For example, the QCD sum rules for the nucleon provided a connection between the 
nucleon mass and the scalar quark condensate (0|gg|0) |3ISI- A more detailed analysis re- 
quires inclusion of the expectation values of the higher dimension, e.g., the gluon condensate 
{Ol^G^iyC^lO) and the four-quark condensate (0|gggg|0). Two other unknowns of nucleon 
QCD sum rules equations are the residue of the nucleon pole Ag and the effective continuum 
threshold Wq. 

Later the QCD sum rules were applied for investigation of modified nucleon parameters 
in the symmetric nuclear matter [HIH]. They were based on the Borel-transformed dispersion 
relation for the function Ilm{(l) describing the propagation of the system with the quantum 
numbers of the nucleon (the proton) in the nuclear matter. Considering nuclear matter as a 



system of A nucleons with momenta pi, one introduces the vector 

P = ^ , (5) 

which is thus p ~ (m, 0) in the rest frame of the matter. The function 11^(5') can be presented 
as Umil) = nm(5'^, V^(p, 5')) with the arbitrary function (p{p,q) being kept constant in the 
dispersion relations in g^. 

The general form of the function Um can thus be presented as 

n™(g) = q,rUliq\s)+p,r^UQ^s) + ml{q\s). (6) 

The in-medium QCD sum rules are the Borel-transformed dispersion relations for the com- 
ponents ^^^(g^ s) {j = q,p, I) 

nj.tf,.) = i/i=!^^(^*v (7) 

The spectrum of the function Ilm{q) is much more complicated than that of the function 
no(g). The choice of the function ip{p, q) is dictated by the attempts to separate the singu- 
larities connected with the nucleon in the matter from those connected with the properties 
of the matter itself. Since the latter manifest themselves as singularities in the variable 
s = {p + q)^, the separation can be done by putting ip{p, q) = {p + qY and by fixing [71 |Hl E] 

(^(p,g) = {p + qf = s = AEIp (8) 

with Eqf being the total nucleon energy at the Fermi surface. 

By using Eq. ((7j) the characteristics of the nucleon in nuclear matter can be expressed 
through the in-medium values of QCD condensate. The possibility of extension of "pole + 
continuum" model [HI^I to the case of finite densities was shown in [7j [TUj. 

The lowest order of OPE of the Ihs of Eq. ((7j) can be presented in terms of the vector and 
scalar condensates [ZI-[ini- Vector condensates f^ = (M|g*7^g*|M) of the quarks with the 
flavor i (|M) denotes the ground state of the matter) are the linear functions of the nucleon 
densities p„ and pp. In the asymmetric matter both SU(2) symmetric and asymmetric 
condensates 

V, = {M\u{0)^^uiO) + diO)^^d{0)\M) = v; + vi 

and 

vlr^ = (M|w(0)7^u(0) - diOh,d{0)\M) = v;-vi 

obtain nonzero values. In the rest frame of the matter t>^ = v'^Sf^o, v^ = v5f^o, vj^^ = v^~^S^o. 
We can present 

v' = (n|g*7og>)p„ + (p|gSg>)Pp- (9) 

The values (A^|g*7oq'*|A^) are just the numbers of the valence quarks in the nucleons {n\u'jou\n) 
{p\d'yod\p) = 1, {p\u'you\p) = {n\d'yod\n) = 2, and thus 




p„ + 2pp = p - - ^ , i;'^ = 2p„ + Pp = p - + ^ . (10) 




Hence, we obtain 

V{p) = VNP, V^-\P,I3) = PV^N^P (11) 

with 

VN = ^, V^N^ = -1 . (12) 

The Ihs of Eq. (|7j) contains the SU(2) symmetric scalar condensate 

Km{p) = {M\u{0)u{0) + d{0)d{0)\M) , (13) 

and SU(2) asymmetric one 

CUp,P) = {M\u{0)u{0) - d{0)d{0)\M) . (14) 

These condensates can be presented as 

Km{p) = /to + f^{p) , 

f^o = 'tm(O) is the vacuum value, 

k{p) = knp + . . . , kat = {N\uu + dd\N) , (15) 

and 

Cm(p, P) = -PiCNp + ■■■), Cn = {p\uu - dd\p) . (16) 

The dots in the rhs of Eqs. ()15j) and (fT?)|l denote the terms, which are nonlinear in p. In 
the gas approximation such terms should be omitted. The SU(2) invariance of vacuum was 
assumed in Eq. (fT^ . The expectation value njy is related to the vrA^ sigma term aT^N, i-e., 

HH 

^N = ^^ (17) 

with rjia^ standing for the current masses of the light quarks, while the value of a^^j^ can 
be extracted from experimental data on low energy vrA^ scattering ^21 ^] ■ However, the 
value Cat should be calculated under certain model assumptions on the quark structure of 
the nucleon. These were the condensates of dimension c? = 3. 

Turning to the condensates of dimension c? = 4, we find for the gluon condensate 

grr^ip) = {M\^G\0)\M) = go + gip), 
vr 

do = 5'm(0) is the vacuum value, 

9{p) = 9NP+ ... (18) 

with the nucleon expectation value 

g^ = {N\^G'\N) ^ -Im, (19) 

TV y 
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obtained in ^^ in a model-independent way. 

Also, the nonlocal condensate {N\q{0)'y^q{x)\N) provides the contributions oi d = 4 and 
those of the higher dimension. The term of the dimension d = 4 is 

{N\r{Oh,D^q\0)\N) = ^g^^-^^mx2 (20) 

with X2 standing for the second moment of the nucleon structure function [^. 

The condensates of dimensions d = 3,4 determine the leading OPE terms of the order 
g^lng^ and Ing^ of the Ihs of Eq. (|7j). The next to leading terms of the order 1/g^ are 
determined by the expectation values of the four-quark operators. The importance of these 
contributions was analyzed in ^H] • In the gas approximation they can be presented in terms 
of the nucleon expectation values. The latter can be obtained in framework of certain models. 

We shall analyze the sum rules in the gas approximation. It is a reasonable starting point, 
since the nonlinear contributions to the most important scalar condensate K,{p) are relatively 
small at the densities of the order of the phenomenological saturation density po = 0.17 fm~'^ 
of the symmetric matter [TOj. The QCD sum rules approach, applied to the description of the 
nucleon self-energies in the symmetric matter [16^ provided the results which are consistent 
with those, obtained by the methods of nuclear physics. 

In the case of symmetric matter {j3 = 0), the leading OPE terms {d = 3, 4) can be either 
calculated or expressed in terms of the observables. In the asymmetric case we need a model 
for the calculation of the expectation value (n, defined by Eq. (J16j) . In most of the quark 
models the nucleon is treated as a system of the valence quarks and the isospin symmetric 
sea of quark-ant iquark pairs. Under this assumption the condensate (^ is determined by 
the valence quarks. In the models with the nonrelativistic valence quarks ^Af = 1- In niore 
realistic relativistic models (^ < 1 due to the relativistic reduction. 

The calculations of the four-quark condensates require model assumptions on the struc- 
ture of the nucleon. The complete set of the four-quark condensates was obtained in [T7] 
by using the perturbative chiral quark model (PCQM). The chiral quark model, originally 
suggested in jTH|, was developed further in [in]^|2I]- In the PCQM the nucleon is treated as 
a system of relativistic valence quarks moving in an effective static field. The valence quarks 
are supplemented by a perturbative cloud of pseudoscalar mesons, in agreement with the 
requirements of the chiral symmetry. This model reflects the main points of the nowadays 
view on the structure of nucleon. The latter is treated as the system of the valence quark and 
the sea-quarks, contained in pions. The chiral properties are respected. In [T2j the SU(2) 
version of PCQM, which includes only the pions, have been used. Thus, in the calculations, 
which include the terms of the order 1/g^ of OPE obtained in framework of PCQM we must 
use the PCQM value (n = 0.54 [19j. 

In the rhs of the sum rules we describe the nucleon by the relativistic in-medium propa- 
gator [221 

Gn' = q^i^-m-^ (21) 



with the total self-energy 

S = g^7^S, + p^7^ ^ + S, . (22) 

in 

We shall use the QCD sum rules for the calculation of the nucleon characteristics 

S^ = ^_^ , m = ^_^ , i:^=m -m , (23) 

identified with the vector self-energy, Dirac effective mass and the effective scalar self-energy 
— see, e.g., j22]- Two other parameters, to be determined in the approach are the in-medium 
shifts of the effective values of the nucleon residue SX^ and of the continuum threshold SW^. 
We present also the result for the single-particle potential energies 

U = S: + S, . (24) 

We trace the dependence of these characteristics on the total density p and on the asymmetry 
parameter /3 — Eqs. (Hj) and (j21). 

Our approach includes only the strong interactions between the nucleons. The electro- 
magnetic and weak interactions are neglected. 

In the asymmetric matter the characteristics obtain different values for the proton and 
neutron. Considering the asymmetry parameter in the interval — 1 < /3 < 1, one can see 
that the values of any parameter P/v for the proton (p) and neutron (n) are connected as 

PM = P,{-P) . (25) 

We shall present all the values for the proton, using Eq. ()25|) to obtain those for the neutron. 

We carry out the calculations in the gas approximation, including only terms linear in 
p in the Ihs of the sum rules. Thus we can neglect the Fermi motion of the nucleons of the 
matter, which manifest themselves in higher order of density, and put 

s = 4m2 (26) 

in Eq. (jH)) . Having in mind future extensions of the approach, we shall keep the dependence 
on s, using Eq. (pUj) for the specific computations. 

In the simplest Hartree approximation (without multiparticle forces) Up is linear in den- 
sity. If one neglects the Fermi motion of the nucleon, as we did, the same is true for the 
scalar self-energy S^. Thus, the parameters Hy and m* (j23|) exhibit nonlinear behavior due 
to nonzero value of Eg (Eg = in the mean field approximation). In similar way, the charac- 
teristics S^ and m* are not linear in p and (3 in our approach due to the relativity large Ihs of 
the sum rules for the structure 11^. This reflects distribution of the baryon charge between 
the pole and continuum. Another reason of the nonlinear effects (although, numerically less 
important) is the nonlinear structure of the sum rules equations. 

The in-medium nucleon QCD sum rules present nucleon interactions with matter in terms 
of the exchange by uncorrelated quark-antiquark pairs. The latter exchanges correspond to 
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the meson exchanges with the same quantum numbers. Thus the contributions of the lowest 
order OPE terms can be viewed to the exchange by the isoscalar and isovector vector and 
scalar mesons, with the point-like vertices of the interactions, having the standard Lorentz 
structure. Inclusion of the nonlocal effects in the quark condensates corresponds to inclusion 
of nonlocal effects in the vertices of the interactions between the nucleons of the matter and 
the mesons. The four-quark condensates, in which two of the quark operators are averaged 
over the vacuum state correspond to the anomalous Lorentz structures of the vertices of 
the nucleon-meson interactions. In terms of hadronic degrees of freedom the four-quarks 
condensates, in which all the quark operators act inside the nucleons correspond to the 
exchanges by the strongly correlated four-quark systems. Such exchanges can be viewed, 
e.g., as those by two mesons, interacting with the nucleon at the same point. This is 
illustrated by Figs. 1 and 2. 

We compare our results for the nucleon self-energies and potential energies with those, 
obtained by the traditional methods of nuclear physics. We find a promising consistency of 
the results. 

We recall in brief the main points of QCD sum rules approach in vacuum and in nuclear 
matter in Sec. II. We present the general equations in Sec. Ill, providing the solutions in 
Sec. IV. We compare the results to those, obtained by the methods of nuclear physics in 
Sec. V. We summarize in Sec. VI. 

I. GENERAL EQUATIONS 

A. Sum rules in vacuum 

To make the paper self-consistent, we recall the main points of the QCD sum approach in 
vacuum [I1I2]- The function Ilo{q^) (often referred to as "polarization operator") is presented 
as 

no(g) = t /rf^xe^(^^)(0|Tj(x)j(0)|0) (27) 



with j being the three-quark local operator (often referred to as "current") with the proton 
quantum numbers. The usual choice is [2] 

j{x) = e'^'' [u''''{x)Cj,u\x)] i,i,d%x) , (28) 

where T denotes a transpose and C is the charge conjugation matrix. The upper indices 
denote the colors. 

The Ihs of Eq. (j3]) is approximated by several lowest terms of OPE, i.e., nQ(g^) k. 
Hq {(f)- The empirical data are used for the spectral function Im T{Q{q^) on the rhs 
of Eq. (j3)). Namely, it is known, that the lowest laying state is the bound state of three 
quarks, which manifests itself as a pole in the (unknown) point k"^ = iv?. Since the next 
singularity is the branching point k"^ = Wpf^ = {m + itIt^Y^ one can present 

Imn^o(^') = A>/(P-m2) + /^(A;2)^(A;2-W^p\) (29) 



with A^ being the residue at the nucleoli pole. Here j = q,I correspond to the two Lorentz 
structures of Eq. Q, ctq = 1, a/ = m. Thus, Eq. (JH) takes the form 



oo 



-fj,OPE(2\ _ ^n'^J I 1 f f^k 



nr""(g^) = -FS + - w^ dk' ■ (30) 

m'^ — q^ 11 J k-^ — q'^ 

Of course, the detailed structure of the spectral density j^ik?^ cannot be resolved in such 
approach. The further approximations are based on the asymptotic behavior 

P{e) = 1 AUr'^^ik') (31) 

at k"^ ^ |g^| with A denoting the discontinuity. The discontinuity is caused by the logarith- 
mic contributions of the perturbative OPE terms. The usual anzats consists in extrapolation 
of Eq. (jHlll to all the values of fc^ , replacing also the physical threshold Wpf^ by the unknown 
effective threshold Wq, i.e.. 



n J k"^ — q^ 2Tri J k^ — q"^ 



w? 



ph 



Thus Eq. (0)) takes the form 

m^ — g^ 2TTi J k"^ — q^ 



n^o'^^^(g^) = 4^ + 7^ / ° ^ dk- . (33) 



w^ 



The Ihs of Eq. ()33|1 contains QCD condensates. The rhs of Eq. (J33|l contains three 
unknown parameters: m, A^ and Wq. The OPE becomes increasingly true while the value 
|g^| increases. The "pole + continuum" model is more accurate at the smaller values of \q^\. 
Thus one can expect Eq. ()33|) to be true in a certain limited interval of the values of \q^\. 
To improve the overlap of OPE and phenomenological descriptions one usually applies the 
Borel transform defined as 

Bf{^') = lim ^^^(-^^Yfi^') - /(M^), 

Q2 = _g2^ M^ = Q^n (34) 

with M called the Borel mass. It is important in the applications to the sum rules that Borel 
transform eliminates the polynomials and emphasizes the contribution of the lowest state in 
rhs of Eq. (jH^ due to relation 

B -^^ = e-"^'/^' . (35) 

m^ — q'' 



The Borel-transforined form of Eq. ()33|) reads 

oo 

flfPE(M') = xle-^yM^aj + ^ f dee-''/^''AUi''''''{e) (36) 

ZTTZ J 



w- 



and is known as QCD sum rules. Actually, there are two sum rules for the structures IIq 
and Ilg of the function no(g) defined by Eq. Q. 

It appeared to be more convenient to work with Eq. (jHUj) multiplied by the numerical 
factor 327r^. The two sum rules for the nucleon in vacuum can be presented in the form [2] 

Ll{M\W^) = Ao(M2), (37) 

L^(M^iyo^) = mAo(M2) (38) 

with 

Ao(M2) = Xle-"'"'^'" . (39) 

Here Ag = 327r^A^, 

oo 

Ll{M\W^) = 327r^('no'^'^^^(M2) - ^ / rffc2e-'='/*^'Ang'^^^(P)Y 

oo 

Ll{M\Wl) = 327r^fno''^^^(M2)-— ( dk^e-'^"'^'" Mi'f''''{e)]. (40) 

V 27ri J J 



wl 



The Ihs of Eqs. flTTj) and (jHHj) J21IS] have been obtained with the account of the condensates 
of dimension rf = 8, i.e., with the account of the terms of the order 1/g^ in OPE of the 
functions Ho 

Lo(M,H/o) - -j^ + ^^^ + -aL 3]^«> (41) 

f/ o On a cib 272 1 a, o , < 

^MMV^) = 2aM^^i-- + — — -ia^ (42) 

with the traditional notations a = — (27r)^(0|gg|0) = —27t^k,q (we assumed the isotopic 
invariance of vacuum {0\uu\0) = {0\dd\0) = (0|gg|0)), b = {27T)'^go, fil = 0.8 GeV^. Here 
Ei are the functions of the ratio Wq/M"^: Ei = EiiWQ/M"^). They are determined by the 
formulas 



Eo{x) = l-e-\ Ei{x) = l-{l + x)e-\ E2(x) = 1 - ( 1 + x + y j e"". (43) 



The factor 

lnz/VA2 
9 



, 2x InMVA^ ^ ^ 

L{M') = ^—^ (44) 



accounts for the anomalous dimension, i.e., the most important corrections of the order as 
enhanced by the "large logarithms". In Eq. (jl^ A = Aqcd = 0.15 GeV, while u = 0.5 GeV 
is the normalization point of the characteristic involved. Note that the two last terms on 
the rhs of Eq. (|4H) originate from the four-quark condensates {0\ur^uur^u\0) and can be 
expressed through the single term ((0|gg|0))^ only in framework of factorization hypothesis 
P]-j2]. Also, the last term on the rhs of Eq. (j42p is the six-quark condensate, evaluated in 
the same approximation. 

The matching of the Ihs and rhs of Eqs. (jHTj) and (jnH|) have been achieved j21 El in the 
domain 

0.8 GeV^ < M^ < 1.4 GeV^ (45) 

providing the values of the vacuum parameters 

A^ = 1.9 GeV^ , W^ = 2.2 GeV^, (46) 

if m = 0.94 GeV. 

B. Sum rules in nuclear matter 

The OPE terms of the polarization operator in nuclear matter 

n^(g) = t fd^xe'^^^^ (M|Tj(x)j(0)|M) (47) 



contains the in-medium values of QCD condensates. Some of these condensates vanish in 
vacuum, obtaining non-zero values only in medium. The other ones just change their values 
compared to the vacuum ones. 

The spectrum of the function Ilm{Q) is much more complicated, than that of the vacuum 
function no(g^). However, [Zj-fnil, the spectrum of the function Ilm{q'^, s) at fixed value of 
s can be described by the "pole -|- continuum" model at least until we include the terms of 
the order p^ in the OPE of Ilm{q'^, s). 

Using Eqs. (|?I|) - (|^. we present the nucleon propagator as 

Or AT = Ztv ■ 5 5 (^48j 

with E„ and m* defined by Eq. ()23p . while the new position of the nucleon pole is 



,2 



and 

Zn 



[s — 171^)11^/171 — T^l + m 
1 + S^/m 

1 

" (l-S,)(l + E,/m) 

10 



*2 



(49) 



(50) 



Thus, we shall present the dispersion relations for the functions Ui^i^q"^, s) {j = q,p, I) in 
the form 

^'^ m- 

with bq = 1, bp = —T,^, bi = m*. The Borel-transformed sum rules take the form 

Ll{M',W^) = AUM'), (52) 

Ll{M',W^) = -KAUM'), (53) 

Li{M',Wl) = m*AUM') (54) 

with 

Am{M^) = X*^e-"''"^/^' . (55) 

Here 

-^m = -^m ■ ^N (56) 

is the effective value of the nucleon residue in nuclear matter. 

Actually, we shall solve the sum rules equations, subtracting the vacuum effects 

L^{M^ Wl W^) = AUM') - Ao(M2) , (57) 

LPiM^,W^) = -S,A™(M2) , (58) 

L'^iM^, Wl, W^) = m*Am{M^) - mAo{M^) (59) 

with U = U^ - L{. The functions Aq, Ll and L^ are defined by Eqs. ((321), (00 and (021), 
while Lg = 0. 

III. CONTRIBUTIONS TO THE LEFT-HAND SIDE OF THE SUM RULES 

We shall use the current, presented by Eq. (PS|) for the calculation of the function nm(g), 
defined by Eq. ()47|) . The possibility of using it in isospin asymmetric system requires addi- 
tional argumentation, which is given in Appendix A. 

We shall include the terms of the order q^ In g^. In q^ and 1/g^ in the Ihs of the sum rules 
(jSTJ). This corresponds to the terms of the order M^, M^ and 1 in the Ihs of the Borel 
transformed equations (jFTfl - fj^^ . We shall include subsequently the contributions of three 
types. The terms £^(M^) stand for the lowest order local condensates. They contribute 
to g^lng^ and Ing^ terms of Eq. (j51|) . These contributions are similar to simple exchanges 
by isovector vector and scalar mesons between the nucleon on the rhs of Eq. (J3T|) and the 
nucleons of the matter. Figs. 1(a), 2 (a). The terms m^(M^) are caused by the nonlocalities 
of the vector condensate. Fig. 1(d), contributing to the terms of the order Ing^ and 1/q^. 
They correspond to the account of the form factors in the vertices between the isovector 
mesons couple to the nucleons. Fig. 2(d). Finally, uj^{M'^) describes the contributions of 
the four-quark condensates. Figs. 1(b), 1(c), being of the order 1/g^. They correspond to 
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the two-meson exchanges (or to exchanges by four-quark mesons, if there are any) and to 
somewhat more comphcated structure of the meson- nucleon vertices, Fig. 2(c). Thus we 
present the Ihs of Eqs. (l3^ - (fS^ as 

LL = C + «'m + < (60) 

and the Ihs of Eqs. (pTTjl -dHHll 

U = e+u^ + w^ (61) 

with l^ = ii^—i-Q, u^ = ui^, uj^ = uI^—Uq, while Eq and Uq are the corresponding contributions 
in the vacuum case. 

The calculation of the function Umiq"^, s) defined by Eq. (jTTj) is based on the presentation 
of the single-quark propagator in medium 

(M|Tg^(x)g^(0)|M) = G,^(x) - ^{M\q\0)^,q\x)\M)^i:^ - ^{M\q\0)q\x)\M)5^p (62) 

with G{x) = (ix^7^)/(27r^x*^) being the free propagator of the quark in the chiral limit. In 
the OPE calculations we neglect the values of the current masses of the light quarks since 
they are small in our scale. Recall that i denotes the light quark flavor. In the lowest orders 
of OPE two of the quarks are described by the free propagators and only one of the quarks 
is presented by the second or the third term of the rhs of Eq. (|62|) . Inclusion of the higher 
order OPE terms will require the two-quark in-medium propagators. 

At X = the matrix elements in the second and third terms on the rhs are just the 
vector and scalar condensates defined by Eqs. (jHI) and (fTH|) . The contribution of the bilocal 
configurations can be expressed in terms of the higher moments and twists of the nucleon 
structure functions 0. 

The bilocal operators on the rhs of Eq. (|U^ are not gauge invariant. The gauge invariant 
expression, achieved by substitution 



q\x) = q'{0) + x^D^q\0) + ^Xo.X(,Do.Dpq\0) + ■ ■ ■ (63) 

with Da standing for the covariant derivatives, provides the infinite set of the local con- 
densates. The expectation values in Eq. (jH^ depend on the variables (px) and x^. In the 
gas approximation we need only the nucleon matrix elements. For the vector structure the 
general form is 

e;{x) = {M\q\0)^,q\x)\M) = ^ <l>'^{{px) , x^) + tx ^m<l>l{{px) , x^) (64) 

with g*(x) defined by Eq. (jU^ . 

The contribution of the vector structure ()64|) to the polarization operator Ilm{q) can be 
presented as 

C(g) = ^ / ^ {x'e- + {X, 9- + e')x) e^(^^). (65) 
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Here we denoted 

a = a^'j'' (66) 

for any 4- vector a^. For the scalar structure 

KiQ) = --J ^ {M\d{0)d{x)\M)e^^^^\ (67) 

vr* J X" 

A. Local condensates of the lowest dimensions 

To include the local vector and scalar quark condensates of the lowest dimensions, we 
must put X = in Eq. (jMjl and in the matrix element {M\d{0)d{x)\M) in Eq. ()67|1 . Thus 

^1(0,0) = {n\q'-foq'\n)pn + {pWloq'\p)Pp, 

while 

{M\d{0)d{0)\M) = {0\dd\0) + pp{p\dd\p) + pu{n\dd\n). 

We must add also the contribution of the gluon condensate. Technically it can be obtained 
by replacing the vacuum gluon condensate by it's in-medium value — Eqs. (|T8|) and (|T9|l . 

The contributions of these local condensates can be expressed as 

£" = /,^(M^W^>^(p) + /J(M^W^^)^(p), (68) 

F = fS{M',Wl)vP{p,(3), (69) 

i' = fi{M\W^)t\p,P), (70) 

with the dependence on p and [3 being contained in the factors 

v'^{p) = 3p, vP{p,(3) = 3p(l-^Y t'{p,(3)=p{KN + CN(3), (71) 

and the function g{p), given by Eqs. (fT^ and (fT!?|) . The other functions are jTHj 



r2 T,.2 ^ 8^' (^ - ^')M^Eo^ - M^E,^ 



mJ 



3 m L^l^ 



mM',w-) = ''"'"'>- 



L'/o 



nM^wl^- '''''"'''' 



mJ 3 ^4/9 

f^{M',Wl) = -An'M^E,.^. (72) 



The notation E^m {k = 0, 1), Eq. (jlHj) . means that the functions depend on the ratio W^/M' 



ml 
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Now we must find the /^-dependence of tlie nucleon self-energies and also of the parameters 
W^ and X^. Note that there is a simple solution of Eqs. (|^ and (j3i|) 

S,(p,/3) = K{p,0)(l-(^) , (73) 

m*{p,p;KN,CN) = m*{p,0;KN + PCn,0) , (74) 

W^ip,/3) = W^ip,0), (75) 

which is true with the same accuracy as the solutions for the symmetric matter {(3 = 0) [TB] . 
Indeed, assuming that W^ does not change with P (f7^ . we find that the function Am(M^) 
in the Ihs of Eq. (|^ should not depend on /3 ^. This leads to Eqs. (fTIj) and (jTSj). 

B. Inclusion of the nonlocal condensates 

Now we present the functions in the rhs of Eq. ()64|) as 

K{b) = Pn4>l,a{b) + Pv4>l,a{b) ■ (76) 

Due to SU(2) invariance they can be presented in terms of the proton functions 0^/f,-) = <P^p^a{b) 

*i(b) = Pn0a(fe) + Pp0a(fe) (77) 

with i = u,d, i' ^ i. 

Expansion in powers of x^ corresponds to the expansion of the function 11^(5') in powers 
of q^. To obtain the terms of the order g~^ it is sufficient to include two lowest terms of the 
expansions in powers of x^. One can present 0121] 

1 

0Lw((P^),^') = /ciae— (^-)/:(,)(a,x2) (7^ 

with 



faibM.x^) = vUb)i(^) + o^^rn'CaUa). (79) 



Here ?7a(«) = /a('^)O) is the contribution of the quarks with the flavor i to the asymptotic 
of the nucleon structure function 77(a) = Vai*^) + vii'^)^ normalized by the condition 

1 

dari{a) = 3 (80) 



with the rhs presenting just the number of the valence quarks in the nucleon. Thus, expansion 
of the function 

faiPx) =0i((px),O) 

^The numerical solution, which will be considered in Sec. IV provides indeed W^{p, (3) w W^{p,0) with 
the accuracy of 10%. 
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in powers of (px) is expressed through the moments of the distributions rjKa). The moments 
are well known — at least, those, which are numerically important. Also, the first moments 



1 



C = JCaic^)da 





of the distribution .^^(a) were calculated in 2,5, by the QCD sum rules method: 

C = -0.24, e"* = -0.09, ^ = ^" + ^"' = -0.33. (81) 

The moments of the function rjKa) can be obtained by using the equations of motion 
Da'~i"q^{x) = miq\x). Thus, in the chiral limit Pj 



{vic^) = l{{v>')-\{e)) 



Here we denoted 



{Q = ^(e» • (82) 



{f)=f'daf{a) (83) 

Jo 



for any function /(a). 

Note that the nonlocality of the scalar condensate, i.e., of the last term on the rhs of 
Eq. ()fi2|) does not manifest itself in the terms up to 1/g^. The first derivative in (px), as 
well all the derivatives of the odd order vanish in the chiral limit due to QCD equation 
of motion. The next to leading order of the expansion in powers of x^ vanishes due to 
certain cancellations as well as in the case of vacuum [2j for the particular choice of the 
operator j(x) presented by Eq. (J28|l . We neglect the terms, containing more complicated 
gluon condensates, since they are numerically small |2(ij . 

Using Eqs. (J64|l and (ffTjl we obtain for the vector structure (J65|l 

C(g) = lC{q) + K\q), (84) 

2i \ r d'^x 



-/?p/^(y.:-¥.^)e^(''^)p|. 

The further calculations, employing Eqs. (jTHj) and (f?^ are similar to those, carried out 
for the case of symmetric matter ^B] — see Appendix B. 
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Finally, the higher moments and higher twists of the nucleon structure functions provide 
the contributions m* to the Ihs L^ of the sum rules — Eq. (plj) 






{u'{M')+(3u%JM'))p 



in- 



3L4/9 



m 



in- 



3L4/9m [2 



-m^M^Eomiva) + -m^is - m^){i) 



rn'M^EoUiv'^a) - {v^a)) 



(85) 



uP{M^ 



K,(M2)+«,(M2))p 



- ^-^m^M^E,^{r^a^) + ^^m'M^E,^{i) 



u'mAM^ 



57r- 



3L4/9 



3 (m^Ei^ - (s - m')M^E^^) ((77" - r/^)«) 



Q 97 



^6) 



u\M^ 



. 



Here we denote L = L{M'^), Ei^ 
of ^ is determined by Eq. (jST|) . 



(87) 
E,{Wl/M^), Eom = Eo{Wl/M^), Eq. (JH, the value 



C. Inclusion of the four-quark condensates 



The exchange by two quark- ant iquark pairs between the current (j28j) and the matter is 
described in terms of the four-quark expectation values 

H^^ip) = {M\uT^uuT^u\M) and R^^ (p) = {M\dT^duT^u\M) (88) 

with r^'"^ being the basic 4x4 matrices 

r^ = /, r^^ = 75, r^ = 7,, r^ = 7^.75, r^ = ^(7^7.-7.7/.), (89) 

acting on the Lorentz indices of the quark operators. We did not display the color indices 
in Eq. (jHH|) . keeping in mind that the quark operators are color antisjTiimetric. 

In the gas approximation 






(90) 
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The characteristics h^"^ and r^^ can be presented as 

5 
6 



/i^^ = -(^{0\uT''u\0){N\uT'^u\N) + {0\uT'^u\0){N\uT''u\N)) 



+ {N\{uT''u-uT''u)int\N) , (91) 

2 
3 



r^^ = '^(^{0\dT^d\0){N\uT^u\N) + {0\uT^u\0){N\dT^d\N)) 



+ {N\{dT''d-uT^u),.^,\N) . (92) 

The first two (" factorized" ) terms in the rhs of Eqs. ()91|) and ()92|) describe two quark 
operators, acting on the vacuum state, while the other two operators act inside the nucleon. 
The last terms describe the " internal" action of all four operators inside the nucleon. This 
is shown by the lower index "int" . The coefficients 5/6 and 2/3 on rhs of Eqs. (|91|) and (|92p 
present the weights of the color-antisymmetric states ^Hl • 

The contribution of the four-quark expectation values to the in-medium change of the 
polarization operator can be written as 

(n)4, = (n^)4, - (no)4, = ^\T.^^xYh^ + T.^xYr^ 

1 \X,Y X,Y 

+ § f E f^xYh^"^ + E rxYT^A . (93) 

y \X,Y X,Y J 

The matrices fixY and txy can be obtained by using the general expression for the function 
^mig) presented in |2E] 

I^XY = ^ Tr (7J%r^)757"g7^75 , 

TXY = X '^^ (^"'?^/3r^)757"r^7^5 , q = q^r ■ (94) 

Here ^y = 1 if F^ has a vector or tensor structure, while ^y = — 1 in the scalar, pseudoscalar 
and axial cases. The sign is determined by that of the commutator between matrix F^ and 
the charge conjugation matrix C — Eq. (16). 

The products ^xYhf^ obtain nonzero values if the matrices F^ and F^ have the same 
Lorentz structure. In this case all the structures presented by Eq. (jHHI) contribute to (Jlm)'iq- 
The products txy'^'^^ do not turn to zero only if F^ has a vector or axial structure. In the 
latter case F"^ should be an axial matrix as well. In the former case F"^ can be either Lorentz 
scalar or Lorentz vector. 

We denote hf-^ = hf, nxx = /Wx and rf^ = rf, txx = Tx for the similar Lorentz 
structures X and Y. The scalar and pseudoscalar expectation values are Lorentz scalars. 
Thus, their contributions can be expressed through single parameters a^ , a^^, a*^^. The 
latter is true also for the scalar- vector expectation value r^^ . We obtain 

q q 

fJ^s = -^ , f^Ps = ^ , {'rsv)^, = -2qf, . (95) 
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In the other channels the four-quark condensates have more comphcated structure. In the 
vector and axial channels 

V{A) . jy{A)PfiP,y 



(^*^^^^)«^ = ah!i^^9f.u + bl 



HV ■"'- ■ •"'- 77^2 



(rP') =a;rig„„ + 6;raf^. (96) 



' IJ.U ' m 



In the tensor channel 






with 



tf^u,pT = -^ [PpPp9vT + PuPrQ^ip - PpPrQvp " PuPp9p,T ) • (9^ 



m^ 



Using Eq. (jMj) we obtain 



hV _ V -^ ,v Hpq) > a _ A -^ : ,A PiPQ) 



m 
TMrf={-6a^^^-2b^Aq + 2b:"'' 






m^ 



T uT I ^ , '^P{P(1)\ 



l^Th = bU ~ + =^^ • (99) 

We can write Eq. (j93j) in the form 

(n)„=B4 + B^M4 + B'H^ (100) 

with each of the coefficients B^ being the composition of the contributions of the neutrons 
and protons of the matter B^ and B^, i.e., 

B^ = Bip,, + Blp,. (101) 

Using Eqs. ^^ - ^^ we obtain 

B! = {-\al, + ^ag - al, + 4, - ^fej,) - 2 (Sa,^, + b^, + ?,a^^, + bi 
Bf = -bl + bl + 2bl, - 2{bX, - bi 



r,i r,ili 

Bl ^ -2af ^ M . (102) 
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Here we denoted af = hf^, af^ = hf'^^'^. Recall that the lower index "i" stands for the 
neutron or the proton. In the last equahty we defined rf^ = af^p^/m, neglecting the 
nonlocal contributions to the scalar-vector condensate. However in the scalar channel such 
terms provide the contributions of the order 1/g^ to the rhs of Eq. (jl(J(J|) . They will be 
included below — Appendix C. 

The calculations of the contributions to rhs of Eq. p()2j) require the model assumption on 
the structure of the nucleon. The complete set of the four-quark condensate was obtained 
in J7] by using the perturbative chiral quark model (PCQM). The chiral quark model, 
originally suggested in [TH], was developed further in [in]-[2Il- In the PCQM the nucleon 
is treated as a system of relativistic valence quarks moving in an effective static field. The 
valence quarks are supplemented by a perturbative cloud of the pseudoscalar mesons in 
agreement with the requirements of the chiral symmetry. In [TT] the simplest SU(2) version 
of PCQM, which includes only the pions, have been used. The parameters which enter 
Eq. ()102|) are calculated in this version of PCQM. 

There are three types of contributions to the four-quark condensate in the framework of 
this approach. All four operators can act on the valence quarks. Also, four operators can act 
on the pion. There is also a possibility that two of the operators act on the valence quarks 
while the other two act on the pions. Following ^Tj we speak of the "interference terms" in 
the latter case. 

To obtain the contribution of the pion cloud, we need the expectation values of the 
four-quark operators in pions. The latter have been deduced in ^7\ by using the current 
algebra technique. It was shown in ^B] , that in the case of the symmetric matter there is 
a remarkable cancellation of the pion contributions in the function Ilm{q)- Since the pion 
contents of the asymmetric matter manifests itself in the different intensities of vr^ and 
7r~ fields, while the four-quark expectation values in vr"*" and vr" mesons are the same, the 
cancellation takes place in this case as well. This cancellation takes place in any model 
of the nucleon which treats the pion cloud perturbatively. Thus, the contributions of the 
four-quark condensates come from the terms, determined by the valence quarks only and 
from the interference terms. 

This enables us to present Eq. (jHH) as 

hf =2- ^{0\uT''u\0){Ni\{uT''uUNi) + {Ni\{uT''uuT''u)i\Ni) . (103) 



Here the lower index "v" means that the operators act on the valence quarks only. The 
lower index "1" corresponds to the sum of the term in which all the four operators act on 
the valence quarks and the term in which two of the operators act on the valence quarks while 
the other two act on pions. Of course, the first ("factorized") term in the rhs of Eq. ()103|) 
obtains a nonvanishing value only in the scalar case P"^ = /. 

The expectation values of the operators of different flavors, providing nonvanishing con- 
tributions to the rhs of Eq. (jU^ are the scalar- vector condensate 

rf^ = 2 ■ '^{0\dd\0){K,\u^^um + {N,\{ddu^^u)^ntm (104) 
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and 






{N;\(dV'^duT^u)^m (105) 



with X standing for vector or axial structures. The meaning of the lower index " 1" is the 
same as in Eq. ()103|) . 

Note that among the interference terms contributing to the four-quark condensates, there 
is so-called "vertex interference", in which one of the PCQM vertices of the self-energy of 
the valence quark is replaced by the four-quark operator. Some of such terms contain 
the matrix elements {Q\u^^d\T{~) and (0|(i75u|7r+), contributing to the expectation values 
(A^|u75(i(i75u|iV), being connected with the matrix elements {N\dT-^ duT-^ u\N) of all struc- 
tures F"^ by the Fierz transform. On the other hand, they depend on the values of the quark 
masses, since (0|'U75(i|7r~) = — ^ ^ ^^'^ with M^ and F^ denoting the mass and the decay 
constant of pion. In a somewhat straightforward approach one substitutes the current quark 
masses. Following more sophisticated models of the pions J2H1 one should substitute the 
constituent quark masses, thus obtaining the values, which are negligibly small in our scale. 

Using the complete set of the nucleon four-quark expectation values [Ej, we obtain 

(n)4, = (^I,(/3)| + AlM^^f. + <(/?)-^) jiy.P (106) 

with the coefficients A\g{P) being determined by Eqs. p03|) - ()105|) . while 

a = -{27iy{0\uu\0) . (107) 

We use the value (0|'Um|0) = (-241 MeV)^, corresponding to a=0.55 GeV^, employed in jSj. 
Note that a is just a convenient scale for presentation of the results. It does not reflect the 
chiral properties of 114^. 

The q term results mainly as the sum of the expectation value of the product of the four 
M-quark operators, described by the flrst ("factorized" ) term on the rhs of Eq. p()3|l . and that 
of the product of two u and two d-quark operators in the vector channel — Eq. ()105|) . The 
former contributions depend on /?, while the latter do not. Thus, the coefficient Al^ depends 
on f3 strongly. The p term is determined mostly by the expectation value ()105|) in the vector 
channel, with the protons and neutrons providing the equal contributions. This explains the 
weak dependence of the parameter A^^ on /3. The coefficient A^^ is dominated by the flrst 
term on the rhs of Eq. (jl04p . providing stronger dependence on f3. The calculations give 

Alg = -0.11 -0.21P, A^^ = -0.57 + 0.09/3, A^^ = 1.90 - 0.92/3 . (108) 

In the simplifled model of the pion, which does not include the renormalization of the 
quark masses by the interactions, the value of the coefficient Af^ is somewhat different 

Alg = 0.25-0.22/3 , (109) 

while the values of A^^ and A^^ remain unchanged. 
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The contributions of the four-quark condensates to the Ihs of the Borel transformed sum 
rules (jS7j)-(ISni) can be presented in the same way, as for the symmetric case 



2 

m 



2m 
IV. SOLUTIONS OF THE SUM RULES EQUATIONS 



f!^ = -Sn\ fi^ = -8n' a, fi^ = -S7r'ma. (110) 



Now we present the solutions of the sum rules equations, focusing on the functions 
E^(p, /5) and m*(p, /5). We shall include the terms i^,u^ and u^ (Eq. (jHH)) in the succession 
in Ihs of Eqs. (J57|) - (j59j) . Two lowest order OPE contributions to the vector structures q and 
p are presented in terms of the vector and gluon condensates and of the nucleon structure 
functions. These characteristics are either calculated in a model-independent way, or deter- 
mined in the experiments. The lowest order OPE terms in the scalar channel are expressed in 
terms of isotope-symmetric and isotope-asymmetric scalar condensates k^v = {p\uu + dd\p) 
and (n = {p\uu — dd\p) — Eqs. (fT3j) and (fTBj). Here the situation becomes somewhat more 
complicated. 

The expectation value ktv is related to the vrA^ sigma-term crTrAr by Eq. ()17|). The value 
of cr^AT can be extracted from the data on low-energy vrA^ scattering. The procedure consists 
in subtracting the high-order chirality- violating terms a' from the experimental value S^rAr, 
i.e., cr,rAf = ^Triv — o"'. The value a' ~ 15MeV was obtained in [221 tiy dispersion relation 
technique. However, there are some uncertainties in deducing the value of S^^y from the 
experimental data. The canonical value S^jy = (60 ± 8)MeV [I2j is now challenged by the 
higher values 77 ± 6 MeV [13 . Assuming irtu + m^ = 11 eV ^T], we find that Stttv = 64 MeV 
corresponds to k^ = 8. Additional uncertainties emerge, since the true value of the sum 
rriu + TUd may be somewhat larger. 

There is no experimental data on the expectation value Cat. If the nucleon is treated as 
a system of the valence quarks and the isospin symmetric sea of the quark-antiquark pairs, 
the expectation value (n is determined by the contribution of the valence quarks. Thus, the 
reasonable values are Ctv = 1 for the nonrelativistic models and Cat < 1 in the relativistic 
case. Until we include only the leading OPE terms i^ , we can solve the sum rules equations 
for any values of n^- and (n. However the four-quark condensates are obtained in framework 
of a specific perturbative chiral quark model (PCQM). Within this model aj^N = 45 MeV 
PUJ . leading to kn = 8. Thus, to be self-consistent, we must use this value as the basic one 
in the general equations, which include the contributions u^ . Note also, that the values of 
Sttat extracted from the experimental data are correlated with the assumption on the strange 
quark content ^Af = / '^\uu+dd\ ) P^ - '^^^ values S^Af ~ 77 MeV correspond to t/N ~ 0.35, 
with a large part of the nucleon mass being due to the strange quarks. The smaller values 
of E^Af require much smaller values of yN- In PCQM one finds i/n = 0.08 [20], in agreement 
with the smaller values of S^at- The PCQM value (n = 0.54 can be obtained by using the 
results of [T^ . 
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We shall present most of the numerical results for the values 

Kj^ = 8 , Cn = 0.54 . (Ill) 

Anyway, in the linearized version of the sum rules, the nucleon characteristics will be pre- 
sented as the explicit function of the condensates, e.g., of the parameters kn and (n- 

A. Solution of the general equations 

Here we present the solutions of the general equations ()57|) - fj59|) . which are identical to 
Eqs. (IS^ - (jM|) . Recall that we approximate the in-medium condensates by the functions, 
which are linear both in p and (3. However the solutions S„(p, /?) and m*{p, /3) are not linear. 
One can demonstrate this by presenting Eqs. ()5H|1 and ()54|1 as 

rp jl 

E„ = --fi , m* = -fi , (112) 

Lira -L/m 

with the density and (3 dependence of L^ leading to nonlinear behavior of S^, and m* (even 
if we assume W^ = Wq). The nonlinear dependence of rhs of Eq. ()112|1 on W^ also cause 
the nonlinear contributions to S^, and m*. However, they are numerically less important. 

Now we include the terms &,u^ and u^ in succession in Ihs of Eqs. (|K7jl - (|Kn|l . 

1. The role of the lowest order local condensates 

As we have seen, these are the contributions, which contain the vector condensate v{p), 
gluon condensate g{p) and the scalar condensates Km and Cm — Eqs. flTT|) . (fT^ . (HH), (fTHjl . 
Account of these terms corresponds to the one-meson exchanges between the nucleon under 
consideration and the nucleon of the matter, with the point-like structures of the meson- 
nucleon vertices. Fig. 2(a). The solution can be obtained by using Eqs. flUH|) - (f7m) for the 
functions ij. As we have seen, there is a simple solution, expressed by Eqs. (f73|) - (f73|l . The 
procedure of minimization of the difference between the Ihs and rhs of Eqs. (jFfj) -()59 |) indeed 
prefers the values Wm{p, P) ~ Wm{p,0). Thus Eqs. (fTSj) and (fTIj) appear to be true with 
good accuracy. Hence, 



Si^Hp,/5) = S.(p,0)(^l-^^ 

eW(p,/3) = e.(p,o)(^i + ^^ 



m*^P\p, P, kn, Cn) = m*{p, 0, kn + PCn, 0) , 

m*(") (p, /3, kn, Cn) = m* (p, 0, k^ - PCn^O) ■ (113) 

Thus, in the matter with the excess of the neutrons {(3 > 0), we obtain E^") > T,^\ and 
i^*(n) ^ rn*^^\ For example, using the value S^(p, 0) obtained by sum rules approach in 
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il6j (St,(p,0) = 335 MeV), we find S(") - S^f) = 170 MeV, m*(") - m<P^ = 50MeV for the 
neutron matter (/5 = 1) at p = po- 

The minimization procedure choses Wl{f3 = -1) = 2.50 GeV^, Wl{f3 = 0) = 2.30 GeV^ 
WliP = 1) = 2.05 GeV^, E^,") - S^f) = 170 MeV and m*(") - m*(p) = 50 MeV. Thus, 
Eqs. fjll3p work well indeed. However, this approximation is not sufficient for the description 
of the potential energy U{p,f3) (f^ . providing U > for both symmetric and asymmetric 
cases. 

2. The role of the four- quark condensates 

Now we include the four-quark condensate, i.e., we use Eq. (jUUj) for L^ = P + uj^ , with 
uj^ described by Eqs. flllOj) . Inclusion of these terms mimics several contributions on the 
rhs of the sum rules (fF7 j) - (fHH|) . In the condensates uT^uuT^u, presented by Eq. (jl03|) the 
first term, which has a nonvanishing value only if T^ = I, generates a contribution to the q 
structure due to the anomalous Lorentz structure of the interaction between the scalar field 
and the nucleon, caused by the chiral-odd vacuum condensate (0|gg|0). In similar way the 
first term on the rhs of Eq. ()104|1 describes the contribution of the vector meson exchange 
to the scalar structure of the nucleon propagator. The anomalous Lorentz structures emerge 
if the nucleon-meson vertices are treated beyond the lowest order. These contributions are 
illustrated by Fig. 2(b). The second terms on the rhs of Eqs. ()103p and ()104|) and on the rhs 
of Eq. ()105p describe exchanges by the four-quark strongly correlated system (see Fig. 2(c)). 
The condensates presented by Eq. fjl05p have the same values for the proton and neutron. 
Thus their contributions do not depend on f3. 

The differences between the neutron and proton characteristics in the neutron matter at 
p = Po become S^,") - E^f) = 140 MeV, m*(") - m*(p) = -110 MeV. In the simplified model 
for the pion with the current masses of the constituent quarks, where A'lq(P) is given by Eq. 
(Unni), we obtain S^,") - S^p) = 145 MeV, while m*(") - m^^) = -110 MeV at these values of 
p and /3. 

3. The final results 

These contributions come from the account of the x-dependence of the expectation values 
of the vector operators (M|g*(0)7^g*(a;)|M) with g*(x) defined by Eq. (jU^ . As we have seen, 
the nonlocality of the scalar condensates is not important in our case. The nonlocality is 
included by putting L^ = & + u^ + u^ with u^ defined by Eq. (jH^ . We use the structure 
functions, obtained in [201 for the calculation of the terms u'^ and m^. 

Account of the nonlocality of the vector condensate corresponds to inclusion of the form- 
factor of the vertex of the interaction between the vector meson and the nucleon of the 
matter - Fig. 2(d). Recall that similar contribution for the effective scalar meson exchanges 
vanishes in our approximation - Sec. III.B. 
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We find the dependence of tlie nucleon vector self-energy S^ and of the effective mass 
m* on density of the matter and on the asymmetry parameter (3 and show the results in 
Figs. 3 and 4. For example, the differences between the neutron and proton characteristics 
in the neutron matter at p = po are S^,") - S^f) = 110 MeV, m*(") - m*(p) = -70 MeV. In 
the simplified model for the pion with the current quark masses, where A^ {(3) is given by 
Eq. (Unni) we obtain S^,") - S^f) = 115 MeV, m*('') - m*(p) = -65 MeV at these values of p 
and p. The minor change is due to the small change in the /? dependence of the condensate 

Due to the nonlinear character of Eqs. (|Ff j) - (j3^ the self-energies T,y{p, /3) and S*(p, /5) = 
m* (p, (3) — m could have been nonlinear in both p and (3. The nonlinear behavior of these 
characteristics with p manifests itself explicitly. However, the dependence on (3 appears to 
be linear in framework of the accuracy of our computations (see the next Subsection). Thus 
both Tjy and E* can be approximated by linear functions of (3 

^,{p^(3) = Uv,{p)+(3tMp) 
Po \ 

S:(p,/?) = -(5i(p)+/?r,^2(p)) (114) 

Po V / 

with Tz = 1 for the proton, r^ = — 1 for the neutron. The functions Vi^2(p) and 5'i_2(p) 
are shown in Fig. 5. They can be approximated by polynomials of the second order - see 
Appendix D. 

The single-particle potential energy is expressed by Eq. (j^ . At p = po the neutron- 
proton difference of the potential energy caused by the isovector interaction is AUnp ~ 
38/5 MeV at small f3. In Fig. 6 we show the dependence U{p) for several values of f3 for 
both neutrons and protons. Recall that the potential energy is determined with the lower 
accuracy then the self-energies. 

The nucleon residue A,^ and the spectrum threshold W^ exhibit very weak dependence 
on p. Thus we can assume 

A^(p,/?)^A^(p,0), W^{p,P) = W^{p,0). (115) 

B. Explicit expression for the nucleon parameters in terms of QCD condensates 

We can present an approximate solution of Eqs. (J57p - (j59p . in which the nucleon self- 
energies are expressed in terms of the QCD condensates explicitly. We see from Eq. (jll5|) . 
that W^{p,p) ^ W^{p,0), while Wl{p,0) is close to it's vacuum value W^ Eq. (HFH) t!)]. 
Thus, we can put W^ = Wq in the rhs of Eqs. (|^ - (fS^ and ()112|) . This enables us to present 
the proton characteristics S^ and m* as the explicit functions of the quark condensates and 
of the Borel mass M^ : 



-'v 



l + F<i{p,P) Po 
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TTl 



m 



+ {T^{^^ + f3CM)+T^Ai^{(3) 



Po 



1 



l + F^ip^P) 



(116) 
(117) 



with 



F\p,P) 



T>jv + mT^g^ + T^, + P^, + ^A^^P) 



P_ 
Po 



:ii8) 



We denote T^ = T^iM^) and T] = T^{M^, W^) for other j and introduce for i = q,p,I 



Ti{M\W',) = p,fl{M\Wl 



,'m?/M'^ 



X' 










{k = v,g,n) 



(r = l,2) 



KiM' 



poL 



4q 



^0 



(119) 



with the functions /^ and fig defined by Eqs. (ff^ and pi()|l . It is instructive to present 
the density p in units of the observable saturation density of the symmetric matter po = 
0.17fm-3. 

v,g,K,ul,u2;i = q,p,I) 



Note that the functions T^{M'^) defined by Eq. ([11^ (j 
depend on M^ rather weakly. Thus, approximating 

T.{M^) = C] , (120) 

we can replace the functions THM"^) in the Ihs of Eqs. ()119|1 by the constant coefficients 
Cj. Numerically most important functions T^{M'^) and Tl{M'^) can be approximated by 
the constant values with the errors of about 4% and 7%. The largest errors of about 25% 
emerge in the averaging of the functions T^. This solves the problem of expressing the 
in-medium change of nucleon parameters through the values of the condensates. For the 
proton 

P 1 






Po 



po (1 - ^q) 
1 



•^.) ' 



with 



^, = - (ClvN + mClgN + mCl^ + PmCl^ + mC^^l 



The coefficients in the rhs of Eqs. (J121|) - (jl23j) are 



C^, 



-0.062, CI = 0.011 GeV-\ C^ 



4q 



-0.070, 



Po' 



(121) 
(122) 

(123) 



Cl = -0.074, CI, = 0.008, 

Ci = -0.042 GeV, C^ = -0.042 GeV, C^ = -0.063 GeV 
CP = -0.090 GeV, C^(_) = -0.068 GeV, C^ = -0.095 , 
C^i = 0.094 , C:, = -0.020. 



(124) 
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Note that the dependence of J^q on (3 is very weak. (Recall that the lowest order OPE 
terms in the q structure of polarization operator do not depend on /3.) Thus, in the Ihs of 
Eqs. fll21|) . p22|) only the dependence of the numerators on (3 is important. This explains 
the linear dependence of the self-energies S^ and S* on (3. 

The values of S„, m* and J-'q for the neutron are described by Eqs. ()12H) - ()123p with j3 
changed to —[3. Equations (|121|) . (jl22|) enable to obtain E^ and m* during the successive 
inclusion of condensates of the higher dimensions. If only the leading OPE terms are included 
the values provided by Eqs. ()121|) . ()122p actually coincide with the solutions of Eqs. (jFTj) - 
(jSni)- If all the contributions are included Eqs. ()121|) . p22|) reproduce the values of S„ and m* 
with the accuracy of 15% and 10% correspondingly for the symmetric matter. The precision 
of Eqs. ()12H) and p22|) changes with (3. For the neutron matter these equations provide 
the values S^,") — Y^fi = 120 MeV and m*^") — m*^^^ = —60 MeV, comparing to the values 
of 110 MeV and -70 MeV, obtained in the previous subsection from the general solution of 
Eqs. (|nij)-(Ei. 

V. DISCUSSION 

Now we compare our results to those, obtained by nuclear physics methods. The lowest 
order OPE terms in the Ihs of the sum rules describe mainly the exchanges by the localized 
qq pairs. This corresponds to the vector and (effective) scalar meson exchanges between the 
nucleon and the nucleons of matter. These exchanges have the point-like vertices and the 
standard Lorentz structures. 

Inclusion of the higher order OPE terms corresponds to a more complicated picture of 
the meson exchanges between the nucleons. Turning to the four-quark condensates, we can 
separate the two types of terms. In the "factorized" contributions one of the qq operators 
is averaged over vacuum. In the "internal" terms both qq pairs act inside the nucleons 
, Eqs. dnH), (ES), (Unnil-dM)- The first (factorized) term in the rhs of Eq. dTU!^ describes 
the contribution to the g^7^structures of the polarization operator, which contains the scalar 
expectation value {N\uu\N). In similar way the first (factorized) term in the rhs of Eq. ()104|) 
contributes to the scalar structure of the polarization operator, being proportional to the 
vector expectation value (A^|'U7^m|A^). These terms correspond to the anomalous Lorentz 
structures of the nucleon- meson vertices. As to the "internal" terms, i.e., the last terms on 
the rhs of Eqs. (jHH), (|^. they can be interpreted as the exchanges by two-meson systems 
with their local interactions with the nucleon or as the exchanges by four-quark mesons (if 
there are any). 

Inclusion of the nonlocal vector condensates ^(0)7^g(x) means that the vertices of the 
interactions between the nucleons of the matter and the vector mesons do not have a point- 
like structure, requiring rather description by the formfactors. The nonlocality of the scalar 
condensate does not influence the results in our approach. 

A usual subject of calculation is the difference between the characteristics of neutron and 
proton. If only the lowest OPE terms are included, the vector self-energies are determined by 
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the vector condensates. The neutron-proton difference S^") — T,^\ usually attributed to the 
p meson exchange is 170 MeV at p = po and /3 = 1. Inclusion of the four-quark condensates 
and of the nonlocalities subtract 30 MeV and 28 MeV from this value. The lowest order 
OPE terms provide the difference of the effective masses m*^"^ — m*^^) = 50 MeV at the 
same values of p and f3. Inclusion of the four-quark condensates and of the nonlocalities 
adds (-160 MeV) and 40 MeV. This leads to E^,") - E^f) = 110 MeV and m*^'^) - m<P^ = 
-70 MeV in neutron matter at p = po 

The structure of the equations for S^ amd m* ()121|) . ()122p is similar to that (j^ employed 
in nuclear physics. Recall that in the Hartree approximation dependence of Sj, on the density 
deviates from the linear law due to Sg ()23|) . (The same refers to S* if the nucleon Fermi 
motion is neglected). In our approach the nonlinear behavior of the self energies is due to 
nonzero values of J-'g. 

Now we compare the numerical results. Considering the papers, containing relativistic 
calculations, we can compare the vector and scalar self-energies S^, and E* = m* — m. In 
the case of the works, carried out in the nonrelativistic approximation, we can compare 
the nucleon potential energies ?7"'^. We analyze also the contribution to the parameter, 
conventionally denoted as a^ jSI], which is defined as 

e{po,P) = £(po,0)+/3V + 0(/?^) , (125) 

being thus the lowest order term of (3"^ expansion of the averaged binding energy e per nucleon 
at saturation value of density. 

Of course, we cannot expect very good agreement, since our calculations are carried out 
in the gas approximation. The future and more sophisticated calculations should include the 
scalar and four-quark condensates beyond the gas approximation. This would correspond 
to the account of the renormalization of the nucleon interactions with the matter by the 
particle-hole excitations in rhs of the sum rules. Another reason is that the results should 
be corrected for the effects of antisymmetrization of the total final state wave function 
("exclusion effect") ^. 

The general feature of the relativistic calculations is that they provide the positive value 
of the difference S^,") — S^^^ > in the matter with the neutron excess. Also the proton 
effective masses are above the neutron ones in this case, i.e., m**-"^ — m*^^' < 0. Our 
calculations show the same tendency. As to the quantitative results, our values of S(,") — S^^^ 
and m**-"^ —m*^^' appear to be about twice smaller than those, obtained in [33^. Our value of 
the effective mass splitting is also about two times smaller than the result of [31] but is only 
30% smaller than that of [33]. However, we find somewhat smaller discrepancy with the 
relativistic Brueckner-Hartree-Fock (RBHF) calculations, presented in jSH]. They found 
j^in) _ j](p) ~ 30 MeV and m*(") - m<P'> ^ -25 MeV at /3 = 0.2 (Fig. 3 of ^), while 
our results are E'^'> - S^p) ^ 20 MeV, m*(") - m*(p) ^ -15 MeV. Another RBHF analysis 
P7j provided the results, which are very close to our ones. One can extract the values 
E(,") - S(,P) ^ 80 MeV, and m*{n) - m*{p) ^ -50 MeV at /3 = 0.75 from Figs. 10 and 11 of 
j37j . Our values are 80 MeV and 55 MeV correspondingly. Note, however, that the split of 
the effective masses, obtained in |37| is due to the exchange effects only. 
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The nonrelativistic calculations, carried out in various approaches [^ EH] provide t/*^"-* — 
f/(p) Ri 60MeV at p = po and P = 1. This is consistent with the earlier calculations [JUj . 
Our value is 40 MeV at p = po and /? = 1. 

Another important parameter is the symmetry energy — Eq. ()125|) . Note that we 
calculate the quantity 

As = l ([/Wp„ + f/Wpp) , (126) 

which is the true contribution to the energy per nucleon. Using Eq. ()114|) we can express 
Ae = |(li + Si- /3\V2 + 52)). The value 

A e{po, /5) - A £(po, 0) = f3^Aa^ + 0(/?^) , (127) 

thus being the true contribution of the isovector forces. Our value is Aa^ ^ 10 MeV. This 
is close to the one, obtained in j32j. The /5^ law is true up to /3^ = 1 with 10% accuracy 
in agreement with [HEl IHHl IHnj- To find the total contribution of the potential energy one 
must include the term caused by the exclusion effect, mentioned above. This adds 9 MeV to 
a^ jH21- Including also the contribution of the kinetic energy, we obtain 04 = 29 MeV. The 
various calculations of this parameter provide the values around 30 MeV [HHl^Ill]- Thus, 
our result agrees with those, obtained by nuclear physics methods. 

V. SUMMARY 

We expressed the vector and scalar self-energies of a nucleon in asymmetric nuclear 
matter as function of density p and of the asymmetry parameter /3. We presented the 
nucleon characteristics in terms of the in-medium expectation values of QCD operators. 
The main ingredients are the nonlocal vector condensates (M|'u(x)7o'u(0) ± d{x)'yod{0)\M) , 
the scalar condensates {M\u{0)u{0) ± d{0)d{0)\M) and the four-quark condensates. The 
local vector condensates are calculated easily. The nonlocality of the vector condensates 
is expressed in terms of the nucleon structure functions. The scalar condensate k{p) = 
{M\uu + dd\M) is presented in terms of the observable sigma term. The scalar condensate 
({p,l3) = {M\uu — dd\M) and the four-quark condensates are calculated in framework of 
perturbative chiral quark model (PCQM). 

Although we treat the condensates in the gas approximation, the nucleon characteristics 
are not linear in density. The corresponding equations ()12H) . ()122|) are analogous to the 
equations of nuclear physics beyond the mean field approximation. Also Eqs. (J12ip . (jl22p 
provide explicit expression of the nucleon characteristics in terms of the QCD condensates. 

The successive inclusion of the OPE terms in the Ihs of the sum rules finds direct analogs 
in the meson-exchange description of the interactions of the nucleon in nuclear matter. The 
lowest order OPE terms correspond to the exchanges by the vector and (effective) scalar 
mesons with the point-like vertices of the interactions. The higher order terms correspond to 
the nonlocal structure of nucleon- meson vertices, including the anomalous Lorentz structures, 
and to the exchanges by the strongly correlated four-quark systems. A possible interpretation 
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of the latter contributions is a local two-meson exchange (or the exchanges by the four-quark 
mesons, if there are any jl2]), Figs. 1, 2. 

We obtained the functions S„(p, /5) and m*{p,[3) for all values of /?. We calculated also 
the single-particle potential energy U{p,P). The results are presented in Figs. 3-6. 

Note, that we did not need phenomenological parameters of the nucleon-meson interac- 
tions. We used the condensates which have been either calculated, or expressed in terms of 
the observables. For the four-quark condensates we used, however, the already known input 
parameters of PCQM. While including the condensates of the higher dimension in succession, 
we found the direct analogs of the meson-nucleon exchange mechanisms of nuclear physics. 

Our results for the nucleon self-energies are in reasonable agreement with the results 
of nuclear physics. The value of the symmetry energy is close to the one obtained by the 
nuclear physics methods. 

The work was supported by the grants RFBR - 03-02-17724 and RSGSS - 1124.2003.2. 

APPENDIX A 

Here we discuss the isotopic structure of the nucleon current (f^H|) . 

The current (PH|) does not have a definite isospin. This is not important for the calcula- 
tions in a medium with the isospin 1=0, e.g., in vacuum or in symmetric nuclear matter. In 
these cases the expectation values of uu and dd operators have equal values (the same is true 
for more complicated condensates). This reasoning does not work for asymmetric matter. 

Therefore it is reasonable to start with the pure 1=1/2 current 



j{y) = eabc ■ 



u^{y)Cj,u\y)j,j>^d%y) 



.l=u^{y)C^^d\yh,ru%y) - -^d'^{y)C^y{yHru'^{y) 



(Al) 



However it turns out that even in asymmetric matter an explicit calculation of the polariza- 
tion operator by using the current (Al) provides the same result as that with the original 
loffe current 



This is due to the Fermi statistics of the quarks. In the OPE the current (j^H|l annihilates a 
local three-quark state. The coordinate function of the three quarks is symmetric. The color 
structure is antisymmetric as well as the spin structure of the spin three quark system (clearly 
the current p8|l has spin 1/2). Thus we need the antisymmetric isospin wave function. 
Due to the Fermi statistic of the fermions any symmetric isospin component will vanish 
after the account for all possible permutations of the quarks. This means that only the 
isospin 1/2 component actually survives in the OPE expansion of the polarization operator 
()27p corresponding to the of current ()28p. In other words there is a relation between the 
contributions of the first and two last terms of the rhs of Eq. (Al). This provides a possibility 
to use (28) for the calculations in isospin asymmetric matter. 
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APPENDIX B 



The integrals in the rhs of Eq. (jH^ can be evaluated by using the formula 



x^ 6 






g'2 



^'^,e*(9'-) (51) 



X" 



with q' = q — pa. The last factor is 



^e*^'^'") = -^g'2 ln(-g'2) + ... (52) 

Here the dots denote the terms which will be killed by the Borel transform. 
Thus the Ihs of the sum rules contains the contribution of the form 

X = fda\n{Q'^ + A\a))f{a) 
with Q^ = —q"^, A{a) = q, v^~"^ )~"^ " . ^^Q Qg^^i present 

X = \nQ^ Jdafia) + j da\n ^ ^"^ /(a). (53) 



The first term on the rhs contains the standard logarithmic factor with the cut, running to 
infinity. It is described by the "pole + continuum" model in a standard way. The second 
term contains a finite cut. It describes the singularities in the -u-channel, caused by the 
nonlocal structure of the vector condensate, corresponding to the exchange terms on the rhs 
of the sum rules. We neglect such contributions, thus coming to the Hartree description of 
the nucleon in nuclear matter. 

APPENDIX C 

Here we calculate the contribution of the four-quark condensates to / structure of the 
polarization operator Hm(g). The contribution is provided by the scalar- vector condensate 
r^ defined by Eq. (j92|) with X = S^Y = V . The first two terms determine the "factorized" 
contribution, which can be written as 

4 = -f /^(x,r(x))e^(^^)(0|Jd|0), (CI) 

with the vector ^"(x) defined by Eq. (jMj) . Neglecting for a while the nonlocal contributions, 
we put 6^{x) = ^"(0). In this approximation 

nj = ^,ioc = -|^(0|Jrf|0)(2p, + pO, (C2) 

o y III 
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thus corresponding to Eq. (fTU^ with a^^ = |(0|rfrf|0), af^ = |(0|rfc/|0). The "internal" 
contributions, provided by the second terms in the rhs of Eqs. (|^ and ()104|) are also 
proportional to (pq), with the coefficients {af'^)int determined in ^^l- The factorized local 
terms and the internal terms provide the rhs of the last equality of Eq. (|102j) . 

Taking into account the x-dependence of the vector 6'^{x) in Eq. (CI) we obtain a more 
accurate expression for the "factorized" contribution, which contains the second moments of 
the structure functions and the terms of the x^ expansion of the integrand 



J _ {0\dd\0)m 



-4 



(pq) 



m^ 



{{Vu)Pp + {Vd)Pn) + 2{{ar]u)pp + {aT]d)pn 



1 u^ 1 

+ i{Qpp + {^d)pn) + ^^ — {{Vu)pn + {Vd)pp) 

I m^ q^ 
The notations are explained in Eqs. ^ and (jHH), pi = 0.8 GeV^ (HTl) . 



(C3) 



APPENDIX D 

The functions Vi,2(p) and 5*1^2 (p) introduced by Eq. (|114j) can be approximated by the 
polynomials of the second order: 



P2{x) = bo + bix + b2x'^ 



(Dl) 



with X = pI Pq. The values of the coefficients hi are presented in Table 1. This leads to 
parametrization of the proton potential energy 



U{x) = [x(-0.15 + 0.17x - O.OOx^) + px{-OMx + 0.02x2)]GeV. 



iD2) 



Table The values of the coefficients bi (GeV) of the polynomials P2 (Dl) , which ap- 
proximate the nucleon self-energies. 





bo 


bi 


b2 


Vi 


0.102 


0.015 


0.026 


V2 


-0.036 


-0.008 


-0.011 


Si 


-0.254 


0.150 


-0.114 


S2 


0.034 


-0.035 


0.034 
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1 Figure captions 

Fig.l. Successive inclusion of the OPE terms. The hehx hues denote the current 
SoUd hues denote the quarks. Thick hues show the matter, (a) The local vector and scalar 
condensates are included (contribution of the gluon condensate is not shown), (b) The 
factorized part of the four-quark condensates is included. One of the qq operators acts on 
vacuum, while another qq operator acts inside the nucleon. Small circles stand for averaging 
over vacuum, (c) The internal parts of the four-quark operators are included. All the quark 
operators act inside the nucleon. (d) The nonlocality of the vector condensate is included. 
The dark blobs denote the nonlocal condensates. 

Fig.2. The analogs of the OPE terms in the meson-exchange interaction of the nucleon with 
matter. Figures 2(a-d) correspond to the Figures l(a-d). Solid line denotes the nucleon. 
Thick lines show the matter. Wavy and dashed lines stand for vector and scalar mesons, 
(a) Exchanges by vector and scalar ("effective" ) mesons with the pion vertices, (b) Exchanges 
by these mesons with the anomalous Lorentz structures of the vertices (the dark squares), 
(c) Exchanges by the qq pairs with strong correlations between them, presumably the local 
interactions with two mesons. Double line denotes the two-meson systems, (d) Inclusion 
of the nonlocal structure of the vertices of interaction between the vector mesons and the 
nucleons of the matter. The dark blobs denote the formfactors of the vertices. 

Fig.3. The density dependence of the vector self-energy S^ (solution of Eqs. (|K7jl - (jKnjl ). The 
solid line shows the results for symmetric matter (/? = 0). The dashed and dotted lines show 
the proton and neutron self-energies T,^' and S^,"-* in the neutron matter (P = 1). 

Fig.4. The density dependence of the effective mass m* (solution of Eqs. (jS7j)-(jSni))- The 
solid line shows the results for symmetric matter (/? = 0). The dashed and dotted lines show 
the proton and neutron self-energies m*^^^ and tti**^") in the neutron matter (/3 = 1). 

Fig.5. The density dependence of the functions xVi, XV2, xSi and xS'2 introduced by 
Eqs. flll4j) . X = —. The dotted lines demonstrate the quality of fitting with the simple 
functions on p, as described in Appendix D. 

Fig.6. The density dependence of the single particle potential energy U, Eq. (1241 The solid 
line shows the results for symmetric matter (/? = 0). The dashed and dotted lines show the 
proton and neutron self-energies U^^' and t/*-"-* in the neutron matter (/5 = 1). 
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